Spin-flip scattering and non-ideal interfaces in dirty ferromagnet/superconductor junctions 
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We study the proximity-induced superconducting correlations as well as the local density of states of a ferro- 
magnet, in a ferromagnet/s-wave superconductor heterostructure. We include the effects of spin-flip scattering, 
non-ideal interfaces, and the presence of impurities in the sample. We employ the quasiclassical theory of super- 
conductivity, solving the Usadel equation with emphasis on obtaining transparent analytical results. As our main 
result, we report that in a certain parameter regime the spatial oscillations of the anomalous (superconducting) 
part of the Green's function induced in the ferromagnet by the proximity effect from the s-wave superconduc- 
tor, are damped out due to the presence of spin-flip processes. As a consequence, spin-flip scattering may under 
certain conditions actually enhance the local density of states due to the oscillatory behaviour of the latter in 
ferromagnet/superconductor structures. We also conjecture that the damping could be manifested in the be- 
haviour of the critical temperature (Tc) of the s-wave superconductor in contact with the ferromagnet. More 
specifically, we argue that the non-monotonic decrease of Tc in ferromagnet/s-wave superconductor junctions 
without magnetic impurities is altered to a monotonic, non-oscillatory decrease when the condition 1 > 16rJ^^ 
is fulfilled, where Tsf is the spin-flip relaxation time and h is the exchange field. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.70.Kn 



I. INTRODUCTION 

Proximity structures consisting of ferromagnetic and super- 
conducting materials offer a synthesis between two important 
physical phenomena that may hold the potential for future ap- 
plications in nanotechnology: spin-polarization and dissipa- 
tionless flow of a current. Ferromagnetism is antagonistic to 
conventional superconductors, since the exchange field acts 
as a depairing agent for spin-singlet Cooper pairs. However, 
the proximity effect does not merely suppress the spin-singlet 
superconducting order parameter, but may also induce long- 
ranged spin-triplet correlations under certain circumstances'. 
Much effort has been invested over the last decade to un- 
veil various physical phenomena that occur in ferromag- 
net/superconductor (F/S) heterostructures^. Among the high- 
lights of such phenomena, it is natural to mention the 7r-state 
that is realized in S/F/S structures, which has been studied in- 
tensively both theoretically ^•"'•^ and experimentally*'^. In this 
state, the superconducting order parameters differ in sign in 
contrast to the usual 0-state in S/N/S structure. The transi- 
tion from a 0- to 7r-state may be controlled by the width of 
the feiTomagnet separating the superconductors, thus offer- 
ing a way of manipulating the Josephson supercurrent that 
occurs in such systems. Another way of obtaining a tt- 
state makes use of misaligned exchange fields in S/F het- 
erostructures. This opportunity arises in a variety of sys- 
tems, ranging from superconductors with spiral magnetic 
order**'^'"', thin S/F bilayer a"''^-'^- '^, and so-called feiTomag- 
netic superconductorsi^ii^ where feiTomagnetic and supercon- 
ducting order seem to coexist uniformly. The latter is most 
often interpreted as evidence for triplet pairing in the super- 
conducting sector 

Although various theoretical ideaUzations allow for a rel- 
atively simple approach to F/S heterostructures in the qua- 
siclassical framework, the presence of factors such as non- 
ideal interfaces and both magnetic and non-magnetic impuri- 
ties should be taken into account to obtain more precise agree- 



ment between theory and experiment. A particularly inter- 
esting feature in such hybrid structures is the generation of 
a spin-triplet superconducting component in the ferromagnet 
which survives even in the dirty limit due to a special symme- 
try property which was first suggested by Berezinskii et al. ' 
and later predicted to occur in F/S junctions by'^''^ Bergeret 
et al. This issue has been the subject of intense investiga- 
tions during the past decade (see for instance Refs. 12 and 
references therein). In particular, the role of triplet pairing in 
superconductor/half-metal/superconductor structures has re- 
ceived much attention latel y^"-^''^^-^- ^, much due to the experi- 
mental verification of a Josephson current in such a setup^"*. 

With regard to F/S junctions, two recent publications have 
adressed some aspects of how spin-flip processes affect the 
critical temperature^'' and the density of states-*. Here, we 
will consider two different geometries to study the impact 
of spin-flip scattering and non-ideal interfaces in heterostruc- 
tures involving ferromagnets and superconductors. The ge- 
ometry of the systems we study are given in Fig. [T]ln the 
top figure, we consider a dirty feiTomagnet of width d sand- 
wiched between a ferromagnetic and superconducting reser- 
voir, where the Green's functions are assumed to be in their 
bulk form. In the bottom figure, the feiTomagnetic reservoir is 
replaced with vacuum, effectively leading to a F/S junction. 

In this paper, we study the influence of magnetic impu- 
rities and non-ideal interfaces on the spatial- and energy- 
dependence of the anomalous (superconducting) part of the 
quasiclassical Green's functions induced in the ferromagnet 
by the proximity effect from the s-wave superconductor In 
particular, we investigate how this is manifested in the local 
density of states (LDOS). We present analytical results that 
may elucidate features obtained numerically in Ref. In 
agreement with Ref. [lH we find that spin-flip processes al- 
ters the decay and oscillation length of the anomalous Green's 
function in the ferromagnet. Our main result is that under 
certain conditions, the usual oscillations of the Green's func- 
tion in F/S junctions without magnetic impurities vanish com- 
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pletely when the condition 1 > 16tjj/i is fulfilled, where r^f 
is the spin-flip relaxation time and h is the exchange field. 
Defining F = l/r^i, an equivalent statement is to say that 
the oscillations vanish when the energy associated with spin- 
flip scattering exceeds a critical value Tc = 4/i. As a direct 
consequence, the spin-flip scattering may actually enhance the 
LDOS due to the oscillations of the anomalous Green's func- 
tion in the ferromagnet. A detailed study is performed con- 
cerning how the length scales associated with the decay and 
the oscillations of the Green's function are affected by mag- 
netic impurities. This is important in the context of under- 
standing the behaviour of for instance the Josephson current 
in S/F/S structures, since spin-flip scattering will always be 
present to some degree in real samples. Including such effects 
will presumably yield a more satisfactory quantitative agree- 
ment with experimental data. 



II. THEORETICAL FORMULATION 

A. Quasiclassical theory 

The central quantity in the quasiclassical theory of 
superconductivity is the quasiclassical Green's functions 
.g(pF, R; £, t), which depends on the momentum at Fermi 
level pf, the spatial coordinate R, energy measured from the 
chemical potential e, and time t. A considerable literature 
covers the Keldysh formalism and non-equilibrium Green's 
functions^^'^*^'^''^^. Here we only briefly sketch the theoreti- 
cal structure, for the sake of readability and for establishing 
notation. The quasiclassical Green's functions .g(pF, R; e, t) 
is obtained from the Gor'kov Green's functions G(p, R; e, t) 
by integrating out the dependence on kinetic energy, assuming 
that G is strongly peaked at Fermi level. 
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FIG. 1 : (Color online) The figure shows geometries that will be con- 
sidered in this paper. In the top figure, we consider a F/F/S junction 
consisting of a dirty ferromagnet sandwiched between a ferromag- 
netic and superconducting reservoir where the Green's functions are 
described by their bulk values. In the bottom figure, we consider a 
F/S junction consisting of a dirty ferromagnet connected to a super- 
conducting reservoir. 



ff(PF,R;e,i) = - / d^pG(p,R;e,i). 

TT 



(1) 



This above is typically applicable to superconducting sys- 
tems where the characteristic length scale of the perturba- 
tions present, such as mean-free path and magnetic coherence 
length, is much smaller than the Fermi wavelength. Also, 
the corresponding characteristic energies of such phenomena 
must be much smaller than the Fermi energy eF- The quasi- 
classical Green's functions may be divided into an advanced 
(A), retarded (R), and Keldysh (K) component, each of which 
has a 4 X 4 matrix structure in the combined particle-hole and 
spin space. One has that 



5=1 ^ 



where the elements of g(pF, R; s, t) read 
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(3) 



The quantities g and / are 2 x 2 spin matrices, with the struc- 
ture 



We organize this paper as follows. In Sec. [Ill we estab- 
lish the theoretical framework we will use to treat the F/S hy- 
brid structure. Namely, we employ the Keldysh formalism 
in the quasiclassical approximation to study the Usadel equa- 
tion with appropriate boundary conditions at the interfaces. In 
Sec. [nil we present our results for the spatial- and energy- 
dependence of the anomalous Green's function in the dirty 
ferromagnet, as well as results for the local density of states, 
both without (for reference) and with spin-flip scattering. We 
also discuss how the decay length and oscillation length scales 
of the Green's function are affected by the spin-flip scattering, 
providing transparent analytical results. In Sec. |lV]and|Vl we 
discuss and summarize the main results of the paper. We will 
use boldface notation for 3-vectors, . T . for 8x8 matrices, . 7 . 
for 4x4 matrices, and for 2x2 matrices. 
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Due to internal symmetry relations between these Green's 
functions, all of these quantities are not independent. In par- 
ticular, the tilde-operation is defined as 



/(pf, R; e, t) = /(-Pf, R; ty 



(5) 



The quasiclassical Green's functions 5 (pf, R; £, t) may be de- 
termined by solving the Eilenberger— equation 



ivfV^ = 0, 



(6) 



where E contains the self-energies in the system such as im- 
purity scattering, superconducting order parameter, and ex- 
change fields. The star-product eg) is noncommutative and 
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is defined in Appendix |A] When there is no expUcit time- 
dependence in the problem, the star-product reduces to normal 
multiplication. This is the case we will consider throughout 
the paper The operation p^g inside the commutator should be 
understood p^g = diagj/Sa, pajg. Pauli-matrices in particle- 
hole x spin (Nambu) space are denoted as pi, while Pauli- 
matrices in spin-space are written as Tj. The Green's func- 
tions also satisfy the normalization condition 



3®5 = 1- 



(7) 



The self-energies entering Eq. (|6]l must be solved in a self- 
consistent manner. For instance, a weak-coupling s-wave su- 
perconducting order parameter is obtained by 



d£(/f|(pF,R;£,i))j 



(8) 



where lu^ is the cut-off energy, which may be eliminated in 
favor of the transition temperature. The notation (...) is to 
be understood as an angular averaging over the Fermi surface. 
Once g(pF, R; s, t) has been determined, physical quantities 
of interest may be calculated, such as the electrical current 



N-pev-p 



deTr{(/53eFg^)pJ, 



(9) 



where N-p is the density of states (DOS) per spin at Fermi 
level. Eq. (|9]l also includes the contribution to charge trans- 
port for holes, thus including processes such as Andreev re- 
flection. In the special case of an equilibrium situation, one 
may express the Keldysh component in terms of the retarded 
and advanced Green's function by means of the relation 



gK^(5R-5A)tanh(/3e/2), 



(10) 



where /3 = is inverse temperature. In nonequilibrium sit- 
uations, one must derive kinetic equations for nonequilbrium 
distribution functions in order to specify the Keldysh part^"*. 

The above equations suffice to completely describe for in- 
stance a single superconducting structure, but must be supple- 
mented with boundary conditions when treating heterostruc- 
tures such as F/S junctions. These boundary conditions take 
different forms depending on the physical properties of the in- 
terface, and we proceed to describe possible scenarios in this 
respect. Transport across interfaces in heterostructures may 
in general be characterized according to three particular prop- 
erties: /) the transmission of the interface, //) the resistivity 
of the compounds separated by the interface, and whether 
the interface is spin-active or not. Let us clarify the distinction 
between the two first properties. The transmission of the bar- 
rier (assuming for simplicity a single open transport-channel) 
determines whether one is dealing with a point-contact or tun- 
neling contact, which differ in terms of the likelihood of elec- 
tron transport to occur across the interface. In the Blonder- 
Tinkham-Klapwijk-language^^ the point-contact corresponds 
to low values of Z while the tunneling limit is obtained for 
high values of Z. On the other hand, the resistivity of the 
compounds separated by the interface is unrelated to the trans- 
missivity of the interface, and one may have for instance a tun- 
neling contact with electrodes attached to it that have either a 
large or small resistance. 



The third property determines to what degree the interface 
discriminates between incoming quasiparticles with different 
spins. Zaitsev^^ derived boundary conditions for a clean N/S 
interface, while Kuprianov and Lukichev (KL)"*^ worked out 
simplified boundary conditions in the dirty limit, valid for 
atomically sharp interfaces in the tunneling regime with a low 
barrier transparency. Although the KL boundary conditions 
are strictly speaking not valid for high transparency of the 
barrier, they may be used for qualitative predictions in that 
regime under certain conditions-''^. The most compact way 
of writing the boundary conditions for the Green's functions 
for arbitrary interfaces was introduced by Nazarov^*^. In all 
the preceding references, a non-magnetic (spin-inactive) in- 
terface was assumed. The generalized boundary conditions 
for magnetically active interfaces have also been derived'*^. 
Let us make a final remark concerning the treatment of in- 
terfaces in the quasiclassical theory of superconductivity. We 
previously stated that the present theory is valid as long as 
characteristic energies of various self-energies and perturba- 
tions in the system are much smaller than the Fermi energy. 
At first glance, this might seem to be inreconcilable with the 
presence of interfaces, which represent strong perturbations 
varying on atomic length scales, clearly in stark contradiction 
to the regime of validity of quasiclassical theory. However, 
this problem may be overcome by including the interfaces as 
boundary conditions for the Green's functions rather than di- 
rectly in the Eilenberger equation. 

The KL boundary conditions may be applied for a dirty 
junction in the tunneling limit when the transparency of the 
interface is low, in correspondance with our assumption of a 
weak proximity effect. For the retarded part of the Green's 
function, they read 



2d-f{gd,j;g) 



x—d 



-[9.9f{V)\ 
[9,9s\ 



(11) 



where F (V) corresponds to the F/F/S (F/S) case of Fig. [T] 
The parameter 7 models the interfacial transmission proper- 
ties, and is given by 7 = Ri/Rp where i?i is the interface 
resistance per unit area, while Rp is the equivalent in the dirty 
ferromagnet. In this work, 7 holds the status of a variational 
parameter. A low transparency of the interface amounts to the 
regime 7^1. As previously mentioned, qualitative predic- 
tions may still be obtained using the linearized Usadel equa- 
tions for 7 ~ 1, and even the quantative aspects of the exact 
numerical solution may in some cases be very well modelled 
by this approximation^ Under the assumption of a weak 
proximity effect, we will neglect the depletion of the super- 
conducting order parameter near the interface in order to fa- 
cilitate the calculations and for the sake of obtaining analyt- 
ical results. Moreover, we will use the bulk solution of the 
Green's function in the superconductor. This approximation 
is valid when the superconducting region is much less disor- 
dered than the ferromagnetl. 
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B. Green's functions 

We will consider the dirty limit of the Eilenberger equation 
Eq. (|6]l, which leads to the Usadel equation"^'. This will be 
an appropriate starting point for diffusive systems where the 
scattering time due to impurities satisfies Xt <^ 1, where X 
is the energy scale of any other self-energy in the problem. For 
strong ferromagnets where h becomes comparable to ep, the 
stated inequality may strictly speaking not be valid for X — 
h. Hence, we will restrict ourselves to the regime h ^ ep- 
Below, we will mostly concern ourselves with the retarded 
part of ^(pp, R; e, t), since the advanced component may be 
found via the relation 



9 



t 



(12) 



The Keldysh component is calculated by means of Eq. ( fTOl i in 
a situation of local thermal equilibrium. For a non-equilbrium 
situation, the Keldysh component is found by 



9 



fT-Tg' 



(13) 



where is a matrix distribution function to be determined 
from kinetic equations, while Eq. ( fT2] l still holds. Eq. (fT3] l 
follows from the normalization condition of the Green's func- 
tion, and the distribution function may be chosen as diagonal 
without any loss of generality. In general, we may write 



(14) 



where comparison with Eq. (fTOl i shows that in a thermal equi- 
librium one has Ti = tanh(/3e/2), Tt = 0. 

By isotropizing the Green's function due to the assumed 
frequent impurity scattering, it is rendered independent of pp. 
This isotropic (in momentum space) Green's function satisfies 
the Usadel equation in the ferromagnet: 



DV(gVg) + i[ep3 + M- a,t, g] - 0. 



(15) 



Above, the exchange energy h is accounted for by the ma- 
trix M = diag(/iT3, /iTa), assuming a magnetization in the 
z-direction, while the spin-flip self-energy reads 



o-(R; e) = -^Y] aig{R; e) 



(16) 



where Tjf is the spin-flip scattering time. We have defined the 
matrices a,; = diag(Ti, r^^). The diffusion constant is given 
hy D — Upr/3. Although the Usadel equation in general re- 
quires a numerical solution, an analytical approach is permiss- 
able under certain conditions. In the case of a weak proximity 
effect, one may effectively linearize Eq. (fTSl ). This is a valid 
treatment for low transparency interfaces or close to Tc- In 
this case, Eq. ( fTSl l is expanded around the bulk solution. To 
be definite, let us consider the retarded part of Eq. ( fTSl l which 
has the same form, namely 

DVig^'yf) + i[efs + M - a,f, g^] = 0, (17) 



where a^f is obtained from dst simply by letting g g^. 
Omitting the superscript on the Green's function, we may ex- 
pand it around the bulk solution go as g ~ go + /, where 
go = diag(l, -1) and 



/ = 











One may now multiply out the matrix equation Eq. ( fTTI ). only 
keeping the lowest order terms in the anomalous Green's func- 
tions fapCR-jS). For more compact notation, we define the 
quantities 

/^^ ^ /;^(R; e), {a, (3} =T, i and /'(^' ^ fff, i = s,t. 

(19) 

We will proceed to consider the two distinct cases illustrated 
in Fig. [U The Green's functions in the different reservoirs 
read 



gv = 0, .gp 



1 
1/' 



.9s 



Ic IT2S 
lf9S —Ic 



(20) 



where we have defined c = cosh(i?), s = sinh(i?) with 
■d = atanh(| A|/e). Note that we have set the superconducting 
phase equal to zero, thus considering a gauge where the gap is 
a purely real quantity. 



C. Odd-frequency pairing 

Before moving on to the graphical presentation of our re- 
sults, let us comment on the presence of the Sz = Q triplet 
component of the anomalous Green's function in the ferro- 
magnet. It is well-known that even in the absence of spin-flip 
processes (rjf — > 00), the triplet component is generated in the 
ferromagnet due to the presence of the exchange field h. We 
will later investigate how the magnitude of this triplet compo- 
nent is affected by including spin-flip processes. Also, it is of 
interest to investigate the symmetry properties of the singlet 
and triplet component. Since we are considering the isotropic 
part (with respect to momentum) of the Green's function due 
to the angular averaging in the dirty limit, one would naively 
expect that only the singlet component should be present. This 
is because the singlet anomalous Green's function is usually 
taken to be even under inversion of momentum, while the 
triplet components are taken to be odd under inversion of mo- 
mentum. Recall that inversion of momentum amounts to an 
exchange of spatial coordinates for the field operators, since p 
is the Fourier transform of the relative coordinate r = ri — r2. 
However, another possibility exists that permits for the pres- 
ence of triplet correlations in the ferromagnet, namely a sign 
shift under inversion of energy. This type of pairing has been 
dubbed odd-frequency pairing in the literature, interpreting 
energy as a real frequency. Recall that inversion of energy 
is equivalent to an exchange of time coordinates for the field 
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operators, since e is the Fourier transform of the relative time 
coordinate t = ti — ^2- For a detailed discussion of even- and 
odd-frequency pairing, the reader may consult AppendixiBl 

Let us in passing show that the singlet component is even 
in frequency, while the triplet component is odd in frequency, 
by using the definition in Eq. (IBI6I 1. To do this, we must first 
find the advanced Green's function by exploiting Eq. ( fT2b . 
Direct matrix multiplication leads to 



/ -1 







-1 







-m-e) 

1 




(21) 



From this, one infers that f^{e) — ff'{~e) (even-frequency 
pairing) and //^(e) = —f^{—e) (odd-frequency pairing). 
We have included this short paragraph on even- and odd- 
frequency pairing to emphasize that although even-frequency 
triplet correlations are destroyed in the dirty limit due to 
the isotropization stemming from impurity scattering, odd- 
frequency triplet correlations may persist since these do not 
vanish under angular averaging. Also, it is important to under- 
stand that the triplet pairing we are discussing in the present 
paper is then quite different from the triplet pairing in for in- 
stance Sr2Ru04. In the latter case, the triplet pairing is odd 
in momentum and therefore even in frequency'*^. As a conse- 
quence, superconductivity is highly sensitive to impurity scat- 
tering in Sr2Ru04 and only observed in very clean samples. 

In the problem under consideration, the two important en- 
ergies are the exchange energy h and the BCS gap |A|. As- 
sociated with these energies are two typical length scales: the 
correlation length in the ferromagnet S,f = \/D/h and the 
superconducting coherence length = \J D / {2-kTc), where 
the critical temperature in a weak-coupling superconductor is 
given by |A| ~ 1.76Tc. One may think of is the pene- 
tration depth of the superconducting condensate into the dirty 
ferromagnet. In an experimental situation, one usually has 
/i ^ |A| even for relatively weak ferromagnets, such that 
^ Cs- For the quasiclassical treatment to be valid, one 
must then have |A| <^ h <IC e^. For a Fermi energy of 1 eV, 
it would then be reasonable to consider h in the neighborhood 
of 30 meV and | A| around 1 meV. The effect of D and d may 
be accounted for in the single parameter ex — D /(P, named 
the Thouless energy. This is the relevant energy scale for the 
proximity effect in the case of highly transparent interfaces. 
In the following, we will unless specifically stated otherwise 
fix /i/|A| = 30 to operate within the allowed boundaries of 
our approximations. Since one is often interested in investi- 
gating how various physical properties behave as a function of 
the thickness d of the ferromagnetic layer, it is useful to note 
that for d/^s ~ x, one finds 



27r|A| 
1.76x2 ■ 



(22) 



ni. RESULTS 

We now provide the main results of this paper, namely a 
study of how the triplet correlations and LDOS are affected 
by spin-flip scattering in a F/S and F/F/S junction. When in- 
cluding scattering upon magnetic impurities in the sample, Eq. 
(fTSI l yields the differential equations 



DdUf.±f.) + He±h){f,±f,\ 



1 

2^ 



0. 



(/. ± 3/,) = 0, 



(23) 



Note that we have here assumed an isotropic spin-flip disor- 
der, in contrast to the uniaxial disorder considered in Refs.l25i 
I27II28I We comment more on this in Sec. |IVl Spin-flip pro- 
cesses in combination with a spatially homogeneous exchange 
field do not lead to equal-spin correlations in the ferromagnet, 
although the inclusion of a spin-active barrier will generate 
these components^^'^". Therefore, for the present case of a 
non-magnetic interface, we have that = 0. For the Sz = 
triplet and singlet Green's functions, the general solution of 
Eq. ( I23I 1 reads 



1 



fs 



ATsfh 



(ciK_e-«+^ 



046 



+ C3K_e^ 
where we have defined 

q± = 



C2K+e 



- C4K+e 



q- X 



(24) 



4iTrfe - 2 ± ./I - 16r2/i2 /{2Dt, 



1/2 



K± = 1 ± ^/l - 16t2/i2. 



(25) 



The coefficients {ct] will be determined from the boundary 
conditions of the F/S and F/F/S junctions. These are given by 
Eq. (fTTI) . which may be written in terms of the f± functions. 
For the F/S junction, we have 



ii): djdxf± 



= 0, 



±s - c/± 



(26) 



In the F/F/S case, the condition ii) is still valid while /) must 
be replaced with 



i): d-ydxfi 



x=0 



f± 



x=0 



(27) 



Note that it is implicit here that f± ~ f±- The resulting ana- 
lytical expressions or {ci} read as follows: 



C4 = 



-2s{l + X_/X+) 



C3 



y_ - Y+X./X+ 

C2 = i?C4, Cl = L2C2 + L3C3 + L4C4 



i2.s-CiY+)/X+, 



(28) 
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For convenience, we have defined the following quantities: 

A± = l± lK_/(4Tsf/l), B± = l± lK+/(4Trf/l), 

-Q+ 



_ g+(^ + g-) ^ _i' 

^2 — } TTi J^3 — — 



La 



A+{q+-ij) 



, R 



q+-ip 

q- 

q- 



(29) 



in addition to 



: A±[e'>+''{2c + 2-fdq+) + e-«+%(2c - 2-fdq+)], 
- B±[e'^-'^R{2c+2-fdq^) + e-«-''(2c - 2-fdq^)] 
■A±e-'^+'^{L4 + L2R){2c-2-fdq+). (30) 



In the F/S case, ^ 0, while -0 = (-7^)"^ in the F/F/S 
case. The knowledge of these coefficients completely deter- 
mines the spatial- and energy-dependence of the anomalous 
Green's functions everywhere in the dirty ferromagnet. When 
applying the limit Tjf oo, by making use of 



A±^2S±,+ , B±^2S±,^, q± 



ik-i 



(31) 



one regains well-known results for the scenario without spin- 
flip scattering. Also, it is worth noting that the triplet com- 
ponent vanishes for /i = 0, even in the presence of spin-flip 
scattering. Having calculated the Green's functions, we may 
now study the effect of spin-flip scattering on the LDOS. The 
spin-resolved LDOS is given by 




'J Me) 



±1. 



(32) 



Above, iYp.cr is the DOS at Fermi level for spin species a 
(Np = J2a^F,<yy The deviation SN from flie bulk DOS in- 
side the dirty ferromagnet may be defined as 



6N : 



[N„CR;e)-Np^^]/Np 



(33) 



We also define the normalized LDOS as = 
A^cr(R; e)/(2A'F,(j), such that in the absence of a 
proximity effect, A^ = 1. The oscillations of the LDOS in 
a F/S junction was first reported by Buzdin^"^, and lead to 
observable effects such as the nonmonotonic dependence of 
the critical temperature on the length d of a F/S bilayer'*'' and 
the TT-phase structures that occur in F/S hybrid systems^. 

Consider first Fig. |2] for a plot of the correction 6N to 
the LDOS as a function of position in the dirty ferromagnet. 
Although the maximum amplitude of 5N is suppressed with 
increasing spin-flip scattering F = l/rsf, an interesting fea- 
ture is that the magnitude of the correction (|(5A^|) is in some 
regions actually enhanced due to spin-flip scattering. This 
seems to be a result of the oscillatory behaviour of the LDOS. 
In a N/S junction, where there is no oscillatory behaviour of 



the LDOS, spin-flip scattering would simply cause a reduction 
of the correction 6N. Thus, the role of spin-flip scattering in 
a F/S junction is more subtle than in a N/S junction where it 
simply amounts to a suppression of the LDOS. 

It might seem counter-intuitive that increasing spin-flip 
scattering should increase the correction to the LDOS, since 
the anomalous Green's functions should be suppressed for 
large F. We suggest the following resolvement of this phe- 
nomena. It is clear that the LDOS displays an oscillatory be- 
haviour due to the presence of an exchange field, both with 
and without the spin-flip scattering. However, in the presence 
of spin-flip processes, the period of these oscillations is mod- 
ified. From Fig. |2] it is seen that the peak of the correction 
to SN is suppressed with increasing F. But even though this 
peak becomes smaller, the different periods of oscillation al- 
lows |(5A^(ri)| to outgrow \SN{r2)\ at certain distances from 
the interface, even for Fi > This is a subtle feature 
unique for F/S interfaces in the presence of spin-flip processes 
as compared to N/S junctions. 
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FIG. 2: (Color online) Spatial variation of the deviation from the 
LDOS (SN) for a F/S junction in the presence of spin-flip scattering. 
We have chosen ex/l A| — 1 and 7 = 5. The qualitative features are 
the same for the F/F/S junctions for this particular set of parameters. 
It is seen that for a given energy, increasing spin-flip scattering F — 
1 1 Tsf will increase the oscillation length and reduce the amplitude of 
the Green's function. 

It is interesting to investigate the role of spin-flip scattering 
with regard to the decay and oscillating lengths further. Very 
recently, some aspects of this topic were adressed in Ref. |5T|. 
We here examine in detail some features that occur when spin- 
flip scattering is included in a F/S junction, among them the 
vanishing of the characteristic oscillations of the anomalous 
Green's function / in a certain parameter regime. Consider 
first the case without spin-flip scattering, effectively letting 
Tsf 00. From our previous equations, it is clear that if we 
write 



k± = ^/2i{e ±h)/D = 



ih 



(34) 
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then the real quantities fci ± and fc2,± correspond to the oscil- 
lating part and decaying part of f±, respectively. We ignore 
the e-term since we consider the regime h e, and write 
1^1, ± I = |fe,±| = 1/Cdec- One readily obtains 



^dec — ^osc — 



(35) 



In other words, we recover the well-known fact that the os- 
cillating and decaying length scales of the superconducting 
condensate in the absence of spin-flip scattering are equal^. 
Consider now a finite value of Tsf, where we obtain 



q± 



16T^f/l2 



(36) 



We have neglected the energy term, assuming STsf ^ 1. In 
this case, writing q± — qi^± + 1(?2,± means that qi^± and q2,± 
are associated with the decay and oscillating length, respec- 
tively. We may now distinguish between two cases. If the 
inequality 



1 > \&tW 



(37) 



is satisfied, then q± is purely real. In this case, there are no 
oscillations of the anomalous Green's function in the ferro- 
magnet. Since we assumed that erjf <C 1, this means that the 
exchange field should be sufficiently weak for the vanishing 
of oscillations to take place. For instance, given a spin-flip 
energy of r/|A| = 20 one would need /i/|A| < 5 for the 
oscillations to disappear If the spin-flip energy becomes very 
large, then the oscillations would vanish even for moderate 
exchange fields. This prediction should have easily observ- 
able experimental consequences, manifested for instance in 
the behaviour of the critical temperature as a function of junc- 
tion width d, given that the required parameter regime may be 
experimentally realized. We comment further on this later in 
this paper. 

In the case where 1 < IQr^fh?, q± is no longer purely real, 
and oscillations are again present in f±. It is instructive to 
consider a plot of ^osc and ^dec as a function of spin-flip scat- 
tering to see how the oscillation and decay length are affected 
by these processes. This is done in Fig. [31 where the diver- 
gence of the oscillation length is clearly seen at the critical 
value Fc = 4ft,. Note that as F ^ 0, ^osc and ^dec become 
equal, as previously stated. When Eq. ( |37] ) is satisfied, the de- 
cay length is different for the up- and down-spins. This may 
be seen by introducing ^j^^ = l/q± sls defined by Eq. ( [36l ). 
Our results for F < Fc are consistent with Ref. 25, who re- 
ported that increased spin-flip scattering should increase the 
oscillation length and reduce the decay length. Let us con- 
sider how this effect is manifested in the LDOS, a directly 
measurable experimental quantity. In Fig. |3(b)[ we plot the 
spatial correction to the LDOS for A| — 5 for several val- 
ues of the spin-flip energy. As seen, the oscillations vanish 
as F increases, and are completely absent when Eq. ( [37] i is 
satisfied. The effect we predict should thus be measurable 
via scanning-tunneling microscopy (STM) measurements, by 
probing the LDOS. 

We next consider a surface plot of the correction to the 
LDOS in the {x, e)-plane for a set of parameters that should 



correspond to a quite typical F/S junction. We set A| = 30 
and F/|A| = 0.3. We use the analytical results for the F/S 
case, since the difference from a F/F/S junction is small in the 
low transparency regime. As seen in Fig. 21 the LDOS peaks 
at a; = rf with an amplitude that increases with energy. This 
peak vanishes upon increasing the spin-flip scattering param- 
eter r, corresponding to the effect we predict - namely that the 
oscillations seen in the LDOS of a F/S junction vanish above 
a critical value of the spin-flip scattering energy. 



IV. DISCUSSION 

In F/S junctions without magnetic impurities, it is well- 
known that the critical temperature of the superconductor 
exhibits a non-monotonic dependence on the thickness of 
the ferromagnet layer d (for an extensive treatment of this 
topic, see the review by Buzdin^) . The physical reason 
for the damped, oscillatory behaviour of Tc in such sys- 
tems is not completely understood. It seems reasonable 
to attribute this characteristic feature to the oscillatory be- 
haviour of the Green's functions in the ferromagnet, since 
this non-monotonic behaviour of Tc is not observed in N/S 
junctions^^'^'*. However, the transparency of the barrier may 
also play a key role in the manifestion of oscillations in Tc, 
as argued in for instance Ref. 55. In another experiment, a 
purely monotonic decay of Tc as a function of d was observed 
in a Pb/Ni junction""^. Even though the physical picture is 
not fully resolved, it is highly plausible that the oscillations 
of / are intimately linked to the behaviour of Tc- In this pa- 
per, we have shown that the presence of spin-flip scattering 
may significantly change the qualitative behaviour of / in the 
ferromagnet. Specifically, the oscillations vanish at a critical 
value Fc — Ah. Under the assumption that the characteris- 
tic behaviour of / strongly influences the fashion in which Tc 
decays, it is then clear that one should observe a transition 
from a non-monotonic (damped, oscillating) to a pure mono- 
tonic (damped) decay of Tc upon increasing the concentration 
of magnetic impurities in a sample. Our finding offers a new, 
possible explanation for the experiments where a monotonic 
decay of Tc was observed, namely that the concentration of 
magnetic impurities was such that the critical value Fc was 
exceeded. 

In this paper, we have assumed an isotropic spin-flip disor- 
der, in contrast to the uniaxial magnetic disorder considered 
m Refs.'25''27 28. This leads to somewhat different equations 
for the proximity-induced anomalous Green's function in the 
ferromagnet. Since we have considered a strongly uniaxial 
exchange field, it is implicitly assumed that the presence of 
spin-flip scattering in the plane perpendicular to the exchange 
field still allows for the uniaxial field to be accomodated. In 
the case of strong uniaxial anisotropy, the first of the Usadel 
equations Eq. ( 1231 ) is replaced by^** 

Ddlif, ± /,) + 2i(£ ± h){f, ± /,) - -^(/, ± /s) - 0. (38) 

2Tsf 

Note how the factor 3 in the last term in Eq. (|23] ) now is 
replaced by unity. Following the same line of reasoning that 
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led to Eq. (I36I I, we now obtain 



9± 



1 T 4ir.,f 



(39) 



This quantity is always complex, and hence we always ob- 
tain (damped) oscillations and never a complete supression 
of the oscillations. Thus, the model with isotropic scattering 
and strongly uniaxially anisotropic scattering are qualitatively 
different. The model with isotropic scattering Eq. ( |23] l is ex- 
pected to be most relevant for a weak exchange field, while 
the model Eq. (|40] | is expected to be most relevant for strong 
uniaxial anisotropy^**. Clearly, it would also be interesting to 
investigate a model which interpolates between these two lim- 
its, in order to investigate at what maximum anisotropy in the 
scattering a complete suppression of oscillations can occur. 
One may consider this situation crudely by the following Us- 
adel equation 

Ddlif, ± /,) + 2i(£ ± h){f, ± /s) - (/t ± /3/s) = (40) 

2Tsf 

where we have introduced the parameter f3 to account for the 
unixial ((3 = 1) and the isotropic ((3 ~ 3) case. Again, fol- 
lowing the line of reasoning that led to Eq. ( [36l ). we now find 



q± 



'(1 + /3)/2Tv/(/5-1)V4 



16t2/i2 



2Dnf 



(41) 



From this simple analysis, one would tentatively conclude that 
the case of strong uniaxial anisotropy is special, in that it is 
the only case where one cannot possibly obtain {f3 — 1)^/4 — 
Wr^fh^ — for any finite Tsf and h, which is the requirement 
for suppression of oscillations. For all other values of f3, it 
would be possible to obtain a real q and hence complete su- 
pression of oscillations. Clearly, this matter warrants further 
and detailed investigations. 

In recent publications^^-^*^, Crouzy et al. considered the in- 
teresting problem of a Josephson current in a S/F/F7S struc- 
ture with non-collinear magnetizations in the ferromagnetic 
regions. It was shown that the misorientation angle between 
the ferromagnetic layers could be used to progressively shift 
the junction between a 0- and 7r-state. In deriving their results, 
effects such as spin-flip scattering and non-ideal interfaces 
were omitted for simplicity. Our analytical results account 
for both of these effects, and may thus be useful to obtain a 
quantitatively better agreement for the Josephson effect with 
experimental data by including these phenomena in S/F/F7S 
structures. Work in this direction is now in progress^. 



The main result of this paper is that we show analytically 
how the well-known oscillations of the anomalous Green's 
function induced in the ferromagnet by the superconductor in 
F/S structures without magnetic impurities vanish completely 
above a critical value for the energy associated with spin-flip 
scattering, F. More precisely, we find that the oscillations are 
absent when the condition Vc > 4/i is fulfilled, where h is the 
exchange field. We have argued that this might be experimen- 
tally observable through a transition from a non-monotonic 
(damped, oscillating) to a monotonic (damped) decrease of 
the critical temperature of the junction as a function of the 
thickness of the ferromagnet layer As another consequence, 
we find that increasing spin-flip scattering may actually en- 
hance the local density of states (LDOS) under certain condi- 
tions. This is a quite subtle effect that might seem counter- 
intuitive at first glance. However, it may be understood by 
realizing that the anomalous Green's functions display an os- 
cillatory behaviour in the presence of an exchange field. The 
period of these oscillations is modified in the presence of spin- 
flip scattering. This means that although the amplitude of the 
oscillations decreases for increasing spin-flip scattering, the 
correction to the LDOS may in certain spatial intervals ac- 
tually be greater for large spin-flip scattering than for weak 
spin-flip scattering. This finding may be of importance in or- 
der to correctly interpret LDOS-spectra as obtained from e.g. 
scanning tunneUng microscopy measurements. 
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APPENDIX A: DEFINING THE STAR-PRODUCT 

We here define the star-product which enters the Eilenberger 
equation Eq. (|6|. For any two functions A and B, we have 



(Al) 



where the differentiation operators denote derivation with re- 
spect to the variables T and e in the mixed representation. 
Note that if there is no explicit time-dependence in the prob- 
lem, the star-product reduces to regular multiplication. 



V. SUMMARY 

In conclusion, we have investigated the role of spin- 
flip scattering and non-ideal interfaces in dirty ferromag- 
net/superconductor (F/S) junctions. Our analytical results 
may serve as a basis for calculating other physical quantities 
of interest in F/S multilayers, such as the Josephson current. 



APPENDIX B: EVEN- AND ODD-FREQUENCY PAIRING 

Consider the anomalous Green's function with time-ordering 
operator T, 

/a/3(ri,r2;ti,t2) = -iT{{ija{ri;ti)'ilj0{r2;t2))}, (Bl) 
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which in the mixed representation may be written as 

/a/3(ri,r2;ti,i2) = /a/3(R, r; T, t). 



(B2) 



The Pauli-principle introduces restrictions on this correlation 
function for equal times ti = t2 = t', namely 



/Q/3(ri,r2;i') = -//3a(r2, ri; f')- 



(B3) 



This follows directly from the anticommutation relation for 
the field operators in Eq. ( IB II ). and reflects the fact that the 
fermionic two-particle anomalous Green's function must be 
antisymmetric under exchange of particle coordinates. As- 
sume now for ease of notation that there is no explicit time- 
dependence in the problem and that the system is homoge- 
neous, which allows us to discard the dependence on the 
center-of-mass coordinates R and T. The following argumen- 
tation is valid even if this simplification may not be performed, 
and the equations then hold for each set of points (R, T). By 
a Fourier-transform, we now obtain 

/a/3(p;i) = y"dre-'P7„Mr;<). (B4) 

The Pauli-hmitation Eq. iB3i then reads /q/3(p;0) = 
//3a(— P;0), since equal times give t = 0. This seems to 
indicate that the Green's function must be odd under inver- 
sion of momentum or exchange of spin coordinates. How- 
ever, another possibility exists, as may be seen by Fourier- 
transforming 

/aMP;£) = / dte'''U(3{p;t). (B5) 

In terms of the momentum- and energy-dependent Green's 
functions, the Pauli-principle now dictates that 

J defa/3ip;e) = - J defi3ai-p,e). (B6) 

This gives us two possibilities that are still perfectly com- 
patible with the equal-time restriction: either fai3{P',£) = 
-fl3a{-p;e) or fap{p;e) = -//3a(-p; -e). This is sum- 
marized in the equation 



faf}{p;£) = -//3a(-p; -e), 



(B7) 



which contains all possible symmetry classifications for the 
Green's functions that are compatible with the Pauli-principle. 
These are listed in Tab. |T] Let us also make contact with the 
Matsubara formalism, where the anomalous Green's function 
is defined as 



/a/3(ri,l-2;ri,T2) 



-r{(^„(ri;Ti)^;3(r2;T2))}, (B8) 



and after a Fourier-transformation to the mixed representation 
satisfies 

/^^(p;ic.„)= Tdre^-V^MP;^)' 
Jo 



C(p;r) 



^2 e-'"""^ f'^ 



/a/3(P;iW« 



(B9) 



with T as a complex time, /? as inverse temperature, and fre- 
quencies Un = (2?T,+ l)7r//3. In this technique, one may apply 
the same procedure as for the real-time Green's functions and 
arrive at 



5^[/^^(p;i^„)+/^J-p;ic.„)]=0, 



which also leads to the requirement that 



/^,(p;ic.„) = -/^J-p;-ic.„). 



(BIO) 



(BID 



The real-time retarded and advanced Green's functions may 
be obtained from the Matsubara Green's function by analyti- 
cal continuation as follows ((5^0): 



(B12) 



From Eq. ( IB7I ). one infers that a triplet correlation may be 



TABLE I: Symmetry classifications of the anomalous Green's func- 
tion that are compatible with the Pauli-principle. The "even" and 
"odd" terminology in the notation here points to the symmetry under 
a sign shift in energy, while "singlet" and "triplet" denotes the sym- 
metry under exchange of spins. With these two properties in hand, 
the symmetry classification of momentum is given from the require- 
ment that the entire function must be antisymmetric. 



Spin 


Momentum 


Energy 


Notation 


Odd 


Even 


Even 


Even singlet 


Even 


Odd 


Even 


Even triplet 


Even 


Even 


Odd 


Odd triplet 


Odd 


Odd 


Odd 


Odd singlet 



even under momentum inversion it is odd under energy in- 
version. This scenario has been dubbed odd-frequency pairing 
in the literature. The Pauli-principle can also be expressed by 
the retarded and advanced anomalous Green's functions by 
using Eq. ( IBllb . To see this, we perform an analytical contin- 
uation on the right hand side of Eq. (IBllb . yielding 

lim /^^(p;ic.„) = /^/5(p;£ + i'5) 

= /^^(p;£), (BI3) 
while the same operation on the left-hand side produces 

lim j-/^^„(-p; -i^n)] = -dp; -s iS) 



^fL{~p;-e). (B14) 



Equating the two sides, we finally arrive at 

/^^(p;e) = -/^J-p;-e). 



(BIS) 



Actually, this information is embedded already in the defini- 
tions of the retarded and advanced Green's functions, and Eq. 
( IB15I ) may be verified by direct Fourier-transformation with- 
out going via Eq. ( IBllb . It is also worth underscoring that the 
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Matsubara technique is only valid for equilibrium situations, 
while the Keldysh formalism and the corresponding Green's 
functions is viable also for non-equilibrium situations. The 
distinction between odd- and even-frequency correlations for 
the retarded and advanced Green's functions is now as fol- 
lows: 

Odd-frequency: f^pip; e) = -f^pip; -e), 
Even-frequency: f^pip; e) = f^p{p; -e). (B16) 
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(a) Plot of the characteristic decay and oscillation lengths {^x ) of the anomalous 
Green's function in the presence of spin-flip scattering with A| = 5. 




(b) Spatial variation of the deviation from the LDOS ( SN) for a F/S junction in 
the presence of spin-flip scattering. We have chosen et/I A| = 1, 7 = 5, and 
/i/l A| = 10. The quasiparticle energy has been set to e/| A| = 0.5, but the 
qualitative behaviour is identical for all e < | A | . 



FIG. 3: (Color online) As shown in (b), there are no oscillations 
of the Green's function in the ferromagnet in the parameter range 
r//i > 4 (F = 1/rsf). Inclusion of the energy contribution e brings 
small corrections to this result, but as seen from the figure the con- 
dition Eq. l l37t is a very good approximation. This behaviour is to 
be contrasted with the usual oscillations in F/S junctions as obtained 
without spin-flip scattering. Note that in (a), ^osc formally diverges 
near F = Fc which separates the two parameter regimes where os- 
cillations occur and where they do not occur. 
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FIG. 4: (Color online) Correction to the LDOS (5N) for r//i = 0.01 
(r = l/Tsf). Surface plot of the deviation from the LDOS in the 
(a;, £)-plane for a junction of width d — 2^s and with transparency 
parameter 7 = 5 and exchange field h/\A\ = 30. The most pro- 
truding feature is the peak emerging in the LDOS right at the F/S 
interface, followed by a dip-structure at low energies. 



